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Abstract 

We analyze the influence of noise on magnetic properties of a su- 
perconducting loop which contains three Josephson junctions. This 
circuit is a classical analog of a persistent current (flux) qubit. A 
loop supercurrent induced by external magnetic fleld at the presence 
of thermal fluctuations is calculated. In order to get connection with 
experiment we calculate the impedance of the low-frequency tank cir- 
cuit which is inductively coupled with a loop of interest. We compare 
obtained results with the results in quantum mode - when the three 
junction loop exhibits quantum tunneling of the magnetic flux. We 
demonstrate that the tank-loop impedance in the classical and quan- 
tum modes have different temperature dependence and can be easily 
distinguished experimentally. 

PACS: 85. 25. Am Superconducting device characterization, design, and 
modeling; 

85.25.Hv Superconducting logic elements and memory devices; micro- 
electronic circuits. 

Keywords: Josephson junction; SQUID; Flux qubit; Thermal fluc- 
tuations 
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1 Introduction 



Magnetic flux quantization in superconductors is used, in particular, for re- 
alization of very sensitive magnetometers. One of them is so-called radio- 
frequency (rf) SQUID [Ij. The sensor of rf SQUIDs is a single junction 
interferometer - a Josephson junction which is incorporated in a supercon- 
ducting ring with a sufficiently small inductance L. When an external flux <l>e 
is applied to an interferometer loop, the circulating supercurrent I is induced 
and a flux $i is admitted into the ring: 

= <l>e - IL. (1) 

The phase difference ip across a Josephson junction equals to a normalized 
magnetic flux in the interferometer loop: 

if = 271^ + 27m. (2) 

where $o is the flux quantum, and n is an integer. Since the Josephson 
current is related to the phase difference ip: 

I = IcSimp, (3) 

where Ic is the critical current, Eq. [T]can be rewritten: 

(f = (fe — f^sinip (4) 

where tpe = 27r$e/ $o is the normalized external flux and the constant 

13 = 27f LI J% (5) 

is the normalized inductance of the interferometer. 

From Eq. H] it is clearly seen that the magnetic properties of a single 
junction interferometer are deflned by parameter /3. If /3 < 1 the <fe{f) de- 
pendence is unique (see Fig. 1) and corresponding mode of SQUID operation 
is so-called nonhysteretic. If /? > 1 the (pei'f) dependence is multivalued (see 
Fig. 1) and corresponding mode of SQUID operation is hysteretic. 

A rf SQUID basically consists of a sensor (usually a single junction inter- 
ferometer) inductively coupled to a radio- frequency-biased tank circuit. The 
flux applied to the sensor changes the effective inductance (or/and the 
effective resistance) of the tank-sensor arrangement. Thus, a flux change can 
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be detected as changes in phase (or/and amphtude) of the voltage across the 
tank circuit. 

The classical mode for the single junction interferometer as well as for cor- 
responding rf SQUID in the presence of fluctuations have been investigated 
in detail theoretically as well as experimentally [2-7]. On the other hand the 
quantum mode for this device is difficult to realize. Since an interferometer 
should be hysteretic, it requires the finite LI^ product, and therefore, the 
finite coupling with environment. In order to avoid this problem, a substitu- 
tion of the geometrical inductance to the Josephson one has been proposed. 
Indeed, if the number of Josephson junctions in the loop m > 2 and for 
suitable junctions parameters, a double degenerated state exists at any ge- 
ometrical inductance L. One of the simplest realizations here is a three- 
junction interferometer, which is called also a persistent current (or fiux) 
qubit [8]. Such qubit was fabricated by several teams and quantum regime 
was convincingly demonstrated. 

If external magnetic fiux is 

9? = TT + 2?™. (6) 

the hysteretic interferometer exhibits double-degenerated energy states, see 
Fig. 3. These states correspond to the different directions of the interfer- 
ometer current. If temperature is low enough and for suitable parameters 
of Josephson junctions the magnetic fiux can tunnel between the two poten- 
tial minima. Below we will call the systems under consideration quantum 
if in their dynamics there is quantum tunneling. If, for some reasons, the 
quantum tunneling is suppressed, we will call these systems classical ones. 

It is clear that the presence of the quantum tunneling ensures the absence 
of the hysteresis at the V2e(v^) dependence. On the other hand jumps between 
two energy minima can be originated by the fiuctuations and the hysteresis 
will be washed out. Therefore for both cases considered above the mode of 
the rf SQUID operation will be nonhysteretic. It arises a natural question: 
by analyzing a SQUID output signal is it possible to distinguish between 
quantum mode (interferometer with "quantum leak") and classical mode 
(interferometer in presence of the fiuctuations)? We address this paper to 
that question. 
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2 Classical mode of a flux qubit in the pres- 
ence of fluctuations 



The studied system presents the superconducting circuit (ring) with three 
Josephson junctions, see Fig. 2. We consider the case of small self inductance 
of the ring L — > 0, therefore $e = $j. The phases across each junction in the 
qubit loop ifi satisfy to 

Vl + ^2 + ^■i = Ve- 
in the frame of RSJ model for Josephson junctions the current through 
each junction is: 

hCi dp h d . r r / \ ■ 1 o /r,N 

~2edi^'^'^2eRdi^'^ '^^^^'^ ^^^^''^^ 
where h is the Plank constant, e is the electron charge, Ci and Ri are the 
junctions' capacitances and resistances respectively. We restrict ourself to 
a practical case, when two junctions in the loop are identical: Ci = C2 = 
C, Id = Ic2 = Ic, Ri = R2 = R and the third junction has slightly smaller 
critical current (with the same critical current density) Jcs = ale, 0.5 < a < 1 
and therefore C3 = aC,R^ = R/a. The presence of the 'white noise' is 
given by the fluctuation currents SIi{t) with correlator < 6Ii{t)6Ii{t ) >= 
2kT/Ri6(t — t ) and mean value < 6Ii{t) >= 0. 
In dimensionless units: 

= 2eRIc/fi,uj^t = T (8) 
and for negligible capacitance, Eq. [7] can be rewritten: 

I/Ic = d(pi^2/dT + sinv9i^2 + 5v'i,2 = adips/dr + asinips + 5(P3{t). (9) 

The correlators of Sifi are: 

< 5(^,(r)5(^,(r') >= 2D6,{t - r'), (10) 

where D = kT/Ej, Ej = ^. 

By introducing the phases 9 and x 

29 = (pi + v?2,2x = Vi- V2, 
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and taking into account that ip^, = — '^G, Eqs. [7] can be presented in the 
form: 



dx/dr = — cos6'sinx + 1/2(5(^2 — ^Vi)^ 
(1 + 2a)de/dT = - sin 61 cos X + a sin(v5e - 26) + 6(P3{t) - l/2{6ip2 + Sifi). 

These equations can be reduced to: 

dU 

dx/dr = -— + 6x{t), (11) 

<^^/<^r^-^^% + m (12) 

where ^ 

U {9, x) = ~ cos 9 cos X — -tt cos((/9e — 29) (13) 

is the effective potential and the random forces are: 

< Sx{r)6x{r') >= D6{t - r), < 69{t)69{t') >= ^^-^D(5(r - r'). 

The Langevin equations (11), (12) describe the random motion of the 'par- 
ticle' with coordinates {9, x) in the periodic potential (13), which is a set of 
bistable cells (eight-shaped contours in Fig. 3a). We have numerically inte- 
grated these stochastic equations by Ito's method (see e.g. jS]) for different 
values of parameter a and the strength of the fluctuations D. The typical 
traces of 9{t) and %(''") are shown in Fig. 4. They correspond to random 
motion in bistable potential. Fig. 3. The arrows indicate the switching from 
one unit cell in Fig. 3a to another. With knowledge of 9{t) and xi"^) the 
average circulating current in the ring is obtained as: 

Ii(t>e) = Ic«stn{x + 9)» . (14) 

The averaging << ... >> includes for each value of flux 0e the average 
over time of traces (6'(r), x(t)) and the average over set of 50 traces. 

The calculated in such way the current-fiux curves for different values of 
D and different values of parameter a are presented in Fig. 5 and Fig. 6. 

From Eqs. (11-13) the Fokker-Plank equation for distribution function 
P(x, 9) can be reconstructed: 

9P - 9 .dU Dd- 1 d dU 1 3D 

dr dx dx 2 dx^ l + 2ad9^d9 ' 2 {1 + 2a)^ d9^ ^ ' 
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The Fokker-Plank equation (15) admits the stationary potential solution 
(see [S]) in the special case a = 1, i.e. when all three junctions are identical. 
For a = 1 the analytical solution reads: 

(16) 



N= I I dxdOe-i^^^''') (17) 




Since the potential U is 2tt periodical function of variables x ^"^^ ^ the 
average current in the ring is: 



_ Jq'' Jo'' dxd9sin{x + 0)e p 
C C dxdOe-m-'^^ 



uixfi) 



;i8) 



In Fig. 7 we compare the numerical results (circles) and the /(</)e) obtained 
from the analitical formula (18) (solid line) for the case a = 1 and D = 
0.2. This comparison was used as an additional calibration of our numerical 
procedure, which is working at arbitrary values of a. 



3 The probing of qubit's state in classical and 
quantum modes. 

As we wrote in introduction we probe a qubit with making use of a tank 
circuit by the impedance measurement technique [10]. It was convincingly 
shown [ini [H] that the observable phase difference S{(pe) between tank cur- 
rent Irf and tank voltage Vrf both in classical and quantum modes reads: 

tan5(¥;,)^e-^, (19) 

where is const, which characterizes the inductive coupling of the qubit 
with the tank circuit. By using the results of Section 2 we calculated output 
signal S{(pe) for classical noise affected mode. For different levels of noise D 
(all of those correspond to nonhysteretic regime) the dependencies S{ipe) are 
shown in Fig. 8. 

In the quantum noise free mode (see Appendix A) the nonhysteretic be- 
havior is achieved by the tunneling between two wells. The phase shift 6{(pe) 
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in this case is described by Eqs. (19,22). It is presented in Fig. 8 for the 
same as in the classical case value of the qubit-tank coupling constant 9 and 
experimentally realized qubit's parameters /p$o = 200GHz, A = l.bGHz. 

Comparing classical and quantum modes (Figs. 8,9) we have found that 
in quantum mode the dip on the 5{^Pe) dependence remains constant at wide 
temperature range kT < A. This reflects the fact that tunneling splitting A 
does not depend on the temperature. The depth of the dip is changed with 
temperature - excitations to the upper level depress the value of the average 
qubit's current. For classical mode the situation is rather different. First of 
all for reasonable set of qubit parameters it is impossible to get such nar- 
row and profound dip similar to obtained in quantum mode. Moreover, the 
temperature dependence of 5(v?e) dip demonstrates that its width strongly 
depends on T. Therefore by analyzing the temperature dependence of 5{ipe) 
one can easily distinguish between quantum and classical modes. 

In conclusion we analyzed the temperature dependence of imaginary part 
of impedance for three-junction loop-tank circuit arrangement in quantum 
and classical modes. We argued that impedance for these modes have quite 
different temperature dependencies and, therefore, can be easily distinguished 
experimentally. 

A Quantum mode of a flux qubit 

Since the tunnel splitting in flux qubit is much smaller than the difference 
between upper energy levels, qubits are effectively two-level quantum sys- 
tems. In the two level approximation a flux qubit can be described by the 
pseudo-spin Hamiltonian 



where ax-, <Jz are the Pauli matrices; A is the tunneling amplitude. The qubit 
bias is given hj e = Ip^ofe, where Ip is the magnitude of the qubit persistent 
current and fe = ^e/^o — 1/2. The stationary energy levels can be easily 
found from the Hamiltonian (1201) : 



H{t) = -Aa, - 



■z , 



(20) 



E±{e) = ±Ve^ + A2, 
and the average value of the qubit current at temperature T is: 



(21) 




(22) 
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Note, that the dependence (21) is valid within the narrow interval of $e 
near $e = $o/2, where the potential U (13) is bistable. 
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FIGURE CAPTIONS 
Figure 1. 0(0e) for rf SQUID in nonhystcrctic and hysteretic modes 
Figure 2. The scheme of persistent current (flux) qubit 
Figure 3. a. Contour plot of the potential U{9,x) (13), a = 0.8. b. Bistable 
potential profile along hne 1-2 in Fig3a at 0e = 7r,Q; = 0.8. 
Figure 4. The random motion of phases 9 and x in bistable potential (Fig. 3b) 
for = 7r,a = 0.8 and D = 0.1 

Figure 5. The dependences I{4>e) for a = 0.8 and D = 0.5(1), D = 0.25(2), D — 
0.1(3) 

Figure 6. The dependences /(0e) for -D = 0.1 and a — 0.5(1), a — 0.8(2), a — 
1(3) 

Figure 7. The comparision of numerical (circles) and analitical (solid line) 
calculations. D = 0.2, a = 1 

Figure 8. The phase shift ^(^^e) in classical mode, a = 0.8 and D = 
0.3(1), £) = 0.2(2), L» = 0.1(3) 

Figure 9. The phase shift S{ipe) in quantum mode. T/A = 0.1(1), T/A = 
1(2) 
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Figure 2: The scheme of persistent current (flux) qubit 
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Figure 3: a. Contour plot of the potefgial U{9,x) (13), a — 0.8. b. Bistable 
potential profile along line 1-2 in Fig3a at 0e = 7r,Q; = 0.8. 




Figure 4: The random motion of phas$4 9 and % in bistable potential (Fig. 3b) 
for 0e = 7r,Q; = 0.8 and D = 0.1 



0.9 



0.8 



0.7 



0.6 



0.4 



0.3 



0.2 



0.1 



o»- 

















































3 












m' 

» 




2 \ 
















» : \ 

\ 

\ \ 








• • 




1 


• 






• 








\ 

\ 

» t 

\ : \ 

\ ^ 













0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

(]) /2n 



Figure 5: The dependences /(0e) for a = 0.8 and D = 0.5(1), 
0.25(2), D = 0.1(3) 
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Figure 6: The dependences l[^e) for D = 0.1 and a = 0.5(1), a = 0.8(2) 
1(3) 
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Figure 7: The comparision of numerical (circles) and analitical (solid line) 
calculations. D = 0.2, a = 1 
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Figure 8: The phase shift S{(pe) in classical mode, a — 0.8 and 
0.3(1), = 0.2(2), D = 0.1(3) 
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